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MATHEMATICAL  ANALYSIS  OF  CURRENT  DISTRIBUTION  IN 
CORROSION  CELLS  WITH  CIRCULAR  GEOMETRY 


INTRODUCTION 

In  the  case  of  uniform  corrosion,  existent  local  anodes  and  cathodes  are  distributed 
randomly  ever  the  surface,  and  their  individual  locations  fluctuate  with  time.  In  this 
situation,  the  electrode  has  a “mixed”  potential  [1]  which,  like  the  current  density,  is 
uniform  throughout  the  metal  surface. 

In  many  cases,  however,  anodic  and  cathodic  reactions  occur  on  different  surfaces,  as 
in  the  galvanic  corrosion  of  dissimilar  metals,  or  on  different  parts  of  the  same  surface,  as 
with  localized  geometries  such  as  corrosion  pits  or  crevices.  In  these  instances,  there  is  a 
varying  distribution  of  both  current  and  potential  over  the  metal  surface. 

Waber  and  coworkers  [2  -6]  have  treated  the  current  distribution  in  systems  where 
the  anode  and  cathode  are  coplanar,  parallel  strips  of  infinite  length.  This  method  of 
calculation  has  been  extended  recently  to  two-dimensional  geometries  for  coplanar,  con- 
centric rectangles  [7]  and  circles  [8 * . 

In  the  case  of  circular  geometries,  Gal-Or,  Raz,  and  Yahalom  [8]  have  calculated  the 
total  current  between  coplanar  electrodes  under  a thin  layer  of  electrolyte,  although  these 
authors  did  not  evaluate  the  local  current  along  the  cell  radius.  In  many  instances  of 
localized  corrosion,  the  corroding  metal  has  a circular  geometry.  Examples  include  some 
instances  of  pitting  [9],  corrosion  under  tubercules  on  copper  [10],  and  crevice  corrosion 
under  O-rings  or  washers  [11] . Thus  it  would  be  useful  to  map  the  corrosion  current 
distribution  in  systems  of  circular  geometry. 

The  purpose  of  this  report  is  thus  to  supplement  the  recent  work  of  Gal-Or  and 
coworkers  [8]  by  evaluating  the  distribution  of  current  and  potential  in  systems  of 
circular  electrodes  as  a function  of  variable  cell  dimensions  and  linear  polarization  param- 
eters. Both  bulk  electrolyte  and  thin-electrolyte  layers  will  be  considered. 


MATHEMATICAL  BACKGROUND 

The  mathematical  method  used  by  Waber  and  coworkers  [2  - 6[  is  well  known.  This 
section  will  sketch  the  elements  of  that  method.  The  paper  by  Gal-Or  and  coworkers 
[8]  gives  a detailed  account  of  their  soluvion  for  the  particular  case  of  circular  electrodes 
covered  by  a thin  layer  of  electrolyte.  The  derivation  for  concentric  circular  electrodes 
under  bulk  electrolyte  has  not  been  published,  although  the  derivation  is  quite  similar  to 
the  thin-layer  case.  The  equation  for  the  bulk  case  will  be  derived  in  this  report  for 
completeness.  It  will  be  seen  that  the  current  distribution  for  bulk  electrolyte  is  the 
limiting  case  for  thin  layers. 
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Potential  Distribution 

The  following  will  apply  to  a generalized  geometry  in  which  the  electrode  surfaces 
are  placed  in  the  XY  plane,  as  shown  in  Fig.  1.  No  specific  geometry  is  yet  assumed. 

The  central  assumption  is  that  the  electrostatic  potential  P at  any  point  in  the 
electrolyte  is  given  by  Laplace’s  equation 

b2P  d2P  d2P  _ n 

a*2  + dy2  + a*2  ' (1) 

Laplace’s  equation  is  valid  if  the  following  conditions  are  met: 

(1)  There  are  no  concentration  gradients  in  the  electrolyte  (uniform  composition). 

(2)  The  solution  is  electroneutral. 

(3)  There  are  no  sources  or  sinks  of  ions  in  the  bulk  electrolyte. 

For  more  detail  see  Ref.  12.  To  solve  Eq.  (1),  one  must  solve  a boundary-value  problem 
for  which  the  potential  is  specified  at  the  boundaries  of  the  electrolyte,  including  the 
electrode  surface. 


-ELECTROLYTE 


j-. ANODE  - CATHODE — 


EFFECT  OF 
POLARIZATION 


Fig.  1— Coordinate  system  and  electrode  po- 
tentials along  the  metal  surface  (y  - constant, 
* = 0) 
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Current  Density 

The  current  density  i(A/cm 2)  at  any  point  in  the  electrolyte  is  related  to  the  potential 
by 


bP 

i =-»r,  (2) 

an 

where  a is  the  solution  conductivity  (ohm-1cm-1 ) and  9 P/dn  is  the  normal  gradient  of 
the  potential  (v  cm-1 ) in  the  direction  of  the  interior  of  the  electrolyte.  At  the  electrode 
surface  (z  = 0), 


i(x,  0)  = -o 


dP(x,z) 

bz 


(3) 


z= 0 


Thus  when  Laplace’s  equation  is  solved  for  P(x,  2),  the  corrosion  current  distribution 
follows  from  Eq.  (3). 

Linear  Polarization 


Waber  [2-6],  following  Wagner  [13],  assumed  linear  polarization,  that  is,  that  plots 
of  current  density  vs  electrode  potential  are  linear.  This  assumption,  intended  to  simplify 
the  mathematics,  is  a reasonable  approximation  in  a number  of  cases. 

Figure  1 also  shows  the  potential  distribution  along  the  electrode  surfaces;  Ea  and  Ec 
are  the  potentials  of  the  unpolarized  anode  and  cathode  respectively.  Figure  2 shows  the 
correspond. r>g  linear  polarization  curves.  The  slope  of  the  cathodic  branch  is 


slope  = 


di(x,  0)  i(x,  0)  - 0 


dP{x,  0)  P(x,  0)  - £ 


■: 


J 


1 


1 


J 


Fig.  2— Linear  polarization  curves,  y = constant  (schematic) 
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Rearranging, 


dE 

P(x,  0)  - i(x,  0)  — - = E 
di  c 


(4) 


where  we  have  written  the  simpler  (and  more  conventional)  notation  dE/di  for  dP(x,  0)/ 
di(x,  0),  the  slope  of  the  linear  polarization  curve.  Use  of  Eq.  (3)  for  i(x,  0)  in  Eq.  (4) 
gives 


P(x,  0) 


L 

dE ' 

t 

9 z 

z=0  di. 

(5) 


The  Wagner  polarization  parameter  £ is  defined  [2,  13]  as 


£ = 


(6) 


where  £ has  the  dimensions  of  length.  Use  of  Eq.  (6)  in  Eq.  (5)  gives 

P(x,0)  . =£„. 

A similar  equation  holds  for  the  anodic  branch: 


(7) 


P{x,  0) 


, \ZP(x,z) 


z=0 


= E_ 


(8) 


If  both  branches  of  the  polarization  curve  have  the  same  slope,  £fl  = £c,  = £ , and 


P(x,  0)  -£ 


?>P{x,z) 

- dz  J z=0 


= S(x) , 


where 


(9) 


0 < x < a ; 
a <x < c . 

Figure  3 shows  an  example  of  polarization  data  [14]  which  give  good  linear  plots. 
When  the  current  is  plotted  in  the  usual  logarithmic  way,  the  cathodic  curve  is  seen  to  be 
typical  of  a diffusion-controlled  process,  whereas  the  anodic  curve  is  typical  of  activation- 
controlled  dissolution.  Yet  both  curves  are  linear  in  current  in  the  millampere  range  as 
the  short-circuited  potential  is  approached. 


Tables  1 and  2 list  a variety  of  systems  which  have  been  observed  to  display  linear 
polarization  curves  over  at  least  a portion,  if  not  the  entirety,  of  current  range.  Thus 
the  mathematical  assumption  of  linear  polarization  is  not  always  unrealistic.  More 
serious,  however,  is  the  added  assumption  that  the  anodic  and  cathodic  linear  polariza- 
tion parameters  are  equal. 
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(b)  POTENTIAL  vs  Ag/AgCI  (mV) 

Fig.  3— Potentiostatic  polarization  of  iron  (10-mil  crevice)  and  platinum 
(open)  in  0.06N  NaCl  (a)  Plotted  semilogarithmically  (b)  Plotted  linearly 


Values  of  the  anodic  linear  polarization  parameters  £ a compiled  in  Table  1 range 
from  0 to  5 cm.  Values  of  the  cathodic  parameter  £c  in  Table  2 range  from  approximately 
0.2  to  250  cm.  Comparison  of  £0  and  £c  for  a given  system  shows  that  £c  is  generally 
larger.  Thus  the  assumption  that  £a  = £c  is  generally  not  valid. 

However,  in  the  case  of  a partially  restrictive  third  dimension,  such  as  a crevice  above 
the  anode,  £a  and  £c  are  probably  closer  to  each  other  than  for  totally  open  geometries. 
With  localized  geometries,  as  in  crevices  or  pits,  it  is  well  known  that  the  local  corrodent 
has  a different  composition  than  the  bulk  solution  [17,  18].  Thus  the  conductivity 
a used  to  calculate  £ should  be  that  of  the  internal  electrolyte.  In  the  case  of  iron,  the 
internal  electrolyte  of  a macro  facsimile  of  an  iron  crevice  or  pit  modeled  according  to 
thermodynamic  reasoning  [19]  has  been  found  to  be  approximately  4M  in  FeCl2  [20]. 

For  0.1N  NaCl,  for  instance,  the  anolyte  conductivity  aQ  would  be  approximately  ten 
times  that  of  the  h"!k  catholyte  conductivity  of  [21],  These  considerations  would  tend 
to  bring  anodic  and  cathodic  values  of  £ into  closer  agreement. 
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Thus  the  assumption  that  anodic  and  cathodic  linear  polarization  parameters  are 
equal,  though  not  necessarily  valid,  does  simplify  the  mathematics  to  allow  a first 
estimate  of  current  and  potential  distributions  within  corrosion  cells. 


BULK  ELECTROLYTE 
Mathematical  Analysis 

The  cell  geometry  is  shown  in  Fig.  4.  The  anode  and  cathode  edges  are  coplanar 
circles  of  radii  a and  c respectively.  With  the  circular  geometry,  it  is  convenient  to  recast 
the  Laplace  equation 


d2P  b2P  b2P 

bx2  by2  b z2 

in  cylindrical  coordinates  through  the  usual  relations 

x = r cos  9 

y = r sin  9 

z = z . 


(1) 


(10) 


The  result  is 


b^P  + 1 bP 
br2  + ~ ~dT  + 


(11) 


with  the  term  in  b2P/bO 2 deleted  because  the  potential  is  not  dependent  on  the  angle  0. 
The  general  approach  is  to 


• solve  Eq.  (11)  for  P(r,  z)  subject  to  the  appropriate  boundary  conditions,  which 
wiil  be  listed  below, 


and  then  to 


• evaluate  the  local  current  density  i(r,  0)  at  the  electrode  surface  using  Eq.  (3) 
in  the  form 


i(r,  0)  = 


-o 


bP(r , z) 
bz 


We  next  list  the  appropriate  boundary  conditions: 


02) 


Boundary  condition  1:  No  current  flows  parallel  to  the  metal  surface  across  the 
terminal  edge  of  the  cathode  (r  = c).  Then 


■i 
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■j 
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0 r • a 


Fig.  4 — Oil  geometry  and  electrode  potentials 


dP(r,  z ) 

Hr,  c)  = -a  — - = 0, 

dr 

w J r = r 


dP(ry  z ) 
dr  . 


Boundary  condition  2 No  current  flows  across  ‘no  symmetry  line  r - 0*: 


^iil  - 0 

. Sr  J „o 


•This  boundary  condition  is  required  to  make  P(r,  r)  bounded  at  r = 0 Alternately  we  could  write  the 
boundary  condition  as  lim  P(r,  ^)<  L,  where  L is  finite 
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r 


The  physical  interpretation  is  that  electrons  lost  in  the  anodic  reaction  migrate  to  the 
nearest  cathodic  region  rather  than  taking  the  longer  path  across  the  line  r = 0. 

Boundary  condition  3:  The  potential  P(r,  2)  mast  be  bounded  at  the  upper  physical 
boundary  of  the  electrolyte,  that  is, 

lim  P(r,  z)  < M , (16) 

2 00 


where  M is  some  finite  number.* 

Boundary  conditon  4:  The  final  boundary  condition  describes  the  potential  at  the 
electrode  surfaces  and  is  given  by  Eq.  (8)  in  terms  of  r: 


P(r,  0)  - £ 


3 P(r,z) 


32 


= S(r)  , 


(17) 


J 2=0 


where 


S(r'  ={ 


Ea , 0<r<Q  ; 
E , a<  r < c . 

r ’ 


(18) 


The  solution  to  Laplace’s  equation  in  cylindrical  coordinates,  that  is,  Eq.  (11)  subject  to 
boundary  conditions  1 through  4,  is 


P(r,z)  =M  Ea  4 

ic. 


2 a / a 
xn  c 1 V"  c 


>-xnzlc 


n=1  ‘'o  (*n) 


— 


(19) 


where  Ea  and  Ec  are  the  potentials  of  the  unpolarized  anode  and  cathode,  J0  and  J j 
ere  Bessel  functions  of  order  zero  and  one  respectively,  and  the  xn  are  the  roots  of 
^(x)  which  give  Jl( jc ) = 0.  A detailed  derivation  of  Eq.  (19)  is  given  in  Appendix  A. 
Along  the  electrode  surface,  the  potential  variation  with  radial  distance  is  thus 


""  (!♦!*„) 


YJ°  \n  c 


(20) 


We  next  evaluate  the  current  distribution  along  the  metal  surface  from 

• z)1 


Hr,  0) 


3P(r, ; 
-a  

()Z 


2=0 


(12) 


•Alternately,  there  is  no  current  flow  across  the  upper  boundary  of  the  electrolyte. 


v. 
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Differentiation  of  Eq.  (19)  gives 


or 


The  left-hand  side  of  Eq.  (22)  is  dimensionless  and  may  be  termed  the  “dimensionless 
local  current  density”  i*: 


i* 


Hr,  0) 


(231 


Thus  the  current  distribution  depends  only  on  the  cell  dimensions  a and  c and  on  the 
polarization  parameter  £ or,  more  properly,  on  the  “reduced”  radius  r/c,  the  “reduced” 
cell  dimension  a/c,  and  on  the  “reduced”  polarization  parameter  £/c. 


Numerical  Evaluation 

Figure  5 shows  the  dimensionless  local  current  density  i*  as  a function  of  reduced 
polarization  parameter  £/c  for  a fixed  cell  configuration  of  a/c  = 0.5.  The  points  in 
Fig.  5 were  computed  using  Eq.  (22)  with  n = 500,  except  near  r/c  = 0 and  r/c  npar 
0.5,  where  2000  terms  were  needed  for  sharp  convergence.  The  computer  program  used 
in  the  evaluation  is  given  in  Appendix  B. 

As  seen  in  Fig.  5,  the  local  current  density  is  uniform  across  the  cell  for  the  larger 
values  of  £ , and  there  is  only  a slight  increase  in  corrosion  current  density  at  the 
anode/cathode  junction  for  the  smallest  £ considered. 

The  corresponding  potential  distributions  are  given  in  Fig.  6.  Potentials  were  cal- 
culated from  Eq.  (20)  with  Ec  taken  to  be  0;  n 
6 shows  the  potential  is  uniform  across  the  cell 
and  not  uniform  for  £/c  = 1. 

For  a given  electrolyte  (fixed  conductivity 
value  of  dE/di:  £ = a]dE/di\.  Thus  the  larger 
that  is,  a greate-  change  in  potential  for  a given 


= 100  sufficed  to  evaluate  the  sum.  Figure 
for  £/c  = 100,  nearly  uniform  for  £/c  = 10, 

a),  the  larger  £ corresponds  to  a larger 
£ means  a more  polarizable  electrode, 
current,  as  shown  in  Fig.  7. 


j 

i 


! 


I 


i 

\ 
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Thus  Figs.  5 and  6 show  that  ^he  more  polarizable  the  electrode  system,  the  more 
uniform  the  distribution  of  potential  and  current.  That  is,  the  effect  of  coupling  such 
electrodes  is  to  distribute  the  polarization  effect  throughout  the  corrosion  cell.  For  the 
least  polarizable  case  shown  (£/c  = 1),  anodic  and  cathodic  areas  are  the  least  affected  by 
each  other,  there  is  the  greatest  drop  in  potential  across  the  anode/cathode  juncture,  and 
hence  there  is  some  intensified  attack  at  that  boundary. 

Figure  8 shows  the  effect  of  changing  the  cell  dimensions.  For  a fixed  anode  radius 
and  polarization  parameter,  an  increase  in  cathode  radius  c leads  to  intensified  attack  at 
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Fig.  6 — Potential  distribution  under  bulk  electrolyte  for  different  polariza- 
tion parameters  at  a fixed  cell  size  (o/c  = 0.5).  (Ec  = 0.) 


the  anode/cathode  boundary  for  the  smaller  £ only.  There  is  no  size  effect  for  the  more 
polarizable  electrodes. 


THIN-LAYER  ELECTROLYTE 
Summary  of  Equations 

The  case  of  a thin-layer  electrolyte  has  been  treated  by  Gal-Or  and  coworkers  [8]. 
The  boundary  conditions  are  the  same  as  the  bulk  case,  except  that  the  requirement  that 
the  potential  be  finite  at  the  electrolyte/air  interface,* 

lim  P(r,  z)  < M , (16) 


*or  that  there  be  no  current  across  that  interface. 
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POTENTIAL  E 


Fig.  7 — Effect  of  the  parameter  f on  the  polarization  curve 
(schematic) 


is  replaced  by  the  requirement  that  there  be  no  current  flow  across  the  outer  boundary 
z = b: 


The  solution  to  Laplace’s  equation  subject  to  Eqs.  (14),  (15),  (17)  and  (24)  is: 


E.  McC4FFERTY 


P,r,  0,  = - + 


c2  - a2' 

f-}Bc 


+ 


n=  1 


1 + - x tanh 


IH^) 

~W 


Vxn> 


IT  Jo  {*«  ?) 


(26) 


and’ 


Again  i*  is  dimensionless  and  may  be  termed  the  dimensionless  local  current  density.  It 
depends  only  on  the  cell  dimensions  a/c , b/c , r/c  and  on  the  reduced  polarization  param- 
eter 2/c. 


Both  the  hyperbolic  tangent  and  hyperbolic  cotangent  approach  unity  as  the  argument 
becomes  infinite,  so  that  when  the  thickness  b of  the  electrolyte  layer  approaches  infinity, 
Eqs.  (26)  and  (27)  in  P(r,  0)  and  i*  respectively  revert  back  to  the  expressions  for  the 
bulk  electrolyte,  that  is,  Eqs.  (20)  and  (21). 


Numerical  Evaluation 

Evaluation  of  i*  according  to  Eq.  (29)  is  shown  in  Fig.  9 for  different  electrolyte 
thicknesses  for  a fixed  cell  size  a/c  and  a fixed  reduced  polarization  parameter  2/c.  It  is 
seen  that  the  anodic  current  density  is  concentrated  at  the  anode/cathode  boundary  for 
the  thinner  electrolyte  layers.  The  corresponding  potential  distribution,  calculated  accord- 
ing to  Eq.  (26),  is  given  in  Fig.  10  (with  Ec  again  taken  to  be  zero).  Like  the  current 
distrib’  ‘..an,  the  potential  distribution  is  uniform  only  for  the  thicker  electrolyte  layers. 

The  effect  of  variable  2/c  on  i*  for  a fixed  geometry  (constant  a/c  and  b/c  ) is  shown 
in  Fig.  11.  The  smaller  values  of  2/c  fiesa  polarizable  electrodes)  Wd  to  pronounced  local 
attack  near  the  edge  of  the  anode.  The  corresponding  potential  distributions  are  shown 
in  Fig.  12.  As  might  be  expected,  the  smaller  values  of  2/c  produce  the  sharpest  drops 
in  potential  at  the  anode/cathode  juncture. 

The  effect  of  variable  cell  size  on  i*  is  shown  in  Fig.  13,  where  a,  b.  and  2 are  fixed 
but  c is  allowed  to  vary.  A decrease  in  the  ratio  of  a to  c leads  to  an  increased  local 
attack  near  the  anode  edge. 


The  expression  by  Gal  Or  and  coworkers  for  the  local  current  density  ( Eq  C 3 X > in  Kef  H|  contains 
a misprint.  Also,  our  notation  follows  Waber,  so  that  b and  c are  interchanged  relative  to  Kef  H 
and  £ and  xn  replace  h and  hn  in  Ref  8 
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Fig.  9 — Effect  of  electrolyte  thickness  on  the  current  distribution  for  a 
fixed  ratio  of  anode  radius  to  cathode  radius  (o/c  = 0.5)  and  fixed 
polarization  parameter  (f/c  = 100).  (n  = 500  except  at  r/c  = 0 and 
near  0 5,  where  n - 2000  to  4000.) 


In  summary,  for  a fixed  anode  radius,  local  attack  in  thin  electrolytes  is  promoted  by 
decreasing. 

• The  electrolyte  thickness  b 

• the  polarization  parameter  £ 

or  by  increasing: 

• the  cathode  radius  c. 


1 


1 


j 


Comparison  with  Experiment 

Few  experiments  have  been  done  with  circular  electrodes.  However,  this  current 
mathematical  analysis  is  in  general  agreement  with  the  investigations  of  Rosenfeld  and 
Pavlutskaya  |15,  22].  For  a circular  iron  anode  in  contact  with  a copper  disc  cathode. 
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Fig.  10-Effect  of  electrolyte  thickness  on  the  potential  distribution  for  a 
fixed  ratio  of  anode  radius  to  cathode  radius  (o/c  = 0.5)  and  fixed  polariza- 
tion parameter  (f/c  = 100).  (£c  = 0,  n = 100.) 


there  was  a non-uniform  distribution  of  potential  under  a 165-pm  layer  of  0.17V  NaCl. 

The  difference  in  potential  between  the  points  r = 0 and  r = c was  nearly  300  mV,  with 
most  of  the  potential  drop  at  the  anode/cathode  boundary.  With  bulk  electrolyte,  the 
potential  distribution  was  more  uniform  with  a potential  drop  of  only  50  mV.  Photo- 
micrographs of  the  anode/cathode  boundary  profile  showed  that  penetration  near  the 
contact  was  much  greater  with  the  thinner  electrolyte. 

These  authors  also  reported  current-distribution  curves  for  a Cu/Zn  couple  which 
were  conjoining  rectangles  rather  than  concentric  circles.  However,  for  that  related 
geometry,  the  anodic  attack  was  greatest  near  the  anode  edge  for  the  thinnest  electrolytes 
and  was  uniformly  distributed  for  the  bulk  electrolyte  (0.17V  NaCl). 

In  the  case  of  another  related  geometry,  where  iron  under  crevices  (thin  electrolytes) 
was  in  contact  with  either  iron  or  platinum  under  bulk  electrolyte,  the  attack  was  in- 
tensified near  the  outermost  edge  of  the  iron  anode  [14]. 


TOTAL  CELL  CURRENT 

As  mentioned  earlier,  Gal-Or,  Raz,  and  Yahalom  [8]  have  evaluated  the  total  cell 
current  rather  than  point-by-point  values  along  the  cell  radius.  These  investigators  showed 
that  the  total  current  increases  with: 
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Fig.  11  — Effect  of  polarization  parameter  on  the  current  distribution  for  a fixed  cell 
geometry  ( ale  = 0.5,  blc  - 0 001 ).  (n  = 2000  for  Sic  = 1,  n = 1000  for  Sic  - 10,  and 
n = 500  for  Sic  = 1 00;  in  all  cases,  n = 2000  to  4000  at  0 and  near  ale. ) 


• increases  in  electrolyte  thickness  (up  to  a limit); 

• decreases  in  polarization  parameter; 

• increases  in  anode  radius  (for  a fixed  cathode  radius). 

The  total  anodic  current  was  evaluated  by  integrating  the  local  current  density: 
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b/c  ■ 0 OOl 


0 02  04  06  08 


r c 

Fig.  12  — Effect  of  polarization  parameter  on  the  potential  distribution  fora 
fixed  cell  geometry  ( ale  = 0.5,  bfc  = 0.001 ).  (£*t.  = 0,  ’ * 1 00. ) 


'total, anodic  = J J 0)  r drd6 

Jr=  0 ”tf=0 

Use  of  i(r,  0)  from  Eq.  (27)  for  thin  layers  gives  [8] 


Fig.  13  — Effect  of  cathode  size  on  the  current  distribution  (anode 
size,  solution  thickness,  and  polarization  parameter  are  fixed) 


For  bulk  electrolyte,  the  term  coth  (x„  b/c)  is  replaced  by  1.  An  internal  check  on  the 
current  distribution  can  be  made  by  comparing  the  total  current  determined  by  mathemati- 
cal integration  from  Eq.  (29)  with  the  total  current  determined  by  graphical  integration 
of  the  current  density  distribution  curve  [i*(r/c)  as  a function  of  r/c].  By  definition. 


i*(r,  c) 


\ 

i 


so  that  Eq.  (28)  can  be  written 


Hr,  0) 


2 


-Ec) 


(23) 
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total, anodic 


4?r  - (£  - £ ) 


r-r=a 

L '•&) 


r dr , 


which  can  be  put  in  the  form 


(total, anodic 


4 noc(Ea  -Ec) 


(30) 


(31) 


The  integral  in  Eq.  (31),  labeled  G*,  can  be  determined  graphically  from  current-distribution 
plots  such  as  Figs.  5,  8,  9,  11  and  13.  Comparison  of  Eqs.  (29)  and  (31)  shows  that 


1 total, anodic  _ / £ \ y _ q * 

4 iroc(Ea-Ec)  \c)  “ 


(32) 


where  Z is  the  sum  indicated  in  Eq.  (29). 

Table  3 shows  that  the  total  anodic  current  determined  by  graphical  integration  of 
the  current  density  distribution  curves  is  in  agreement  with  values  determined  by  mathemati- 
cal integration  after  Gal-Or  and  coworkers.  Table  3 also  shows  the  effect  of  the  various 
parameters  on  the  total  anodic  current;  these  trends  have  already  been  reported  by  Gal-Or, 
Ra z,  and  Yahalom  (8]  and  are  listed  at  the  beginning  of  this  section. 

As  a final  check,  the  total  cathodic  current  was  computed  in  some  cases  by  graphical 
integration  of  current  density  distribution  curves  from  r/c  = a/c  to  r/c  = 1.  The  total 
cathodic  current  was  equal  to  the  total  anodic  current,  as  required. 


SUMMARY 

The  distribution  of  current  in  a circular  corrosion  cell  has  been  computed  based  on 
the  assumption  that  the  electrolyte  has  a uniform  composition.  The  mathematical  equations 
are  summarized  in  Taule  4.  Numerical  evaluation  of  the  current-distribution  function 
shows  that: 

• Current  distribution  is  more  uniform  along  the  metal  surface  under  bulk  elec- 
trolytes than  under  thin  layers  of  electrolyte. 

• Under  thin  layers  of  electrolyte,  corrosion  attack  is  intensified  at  the  anode/ 
cathode  junction.  For  a fixed  anode  radius,  the  local  attack  is  promoted  by: 

• a decrease  in  electrolyte  thickness, 

• a decrease  in  the  linear  polarization  parameter, 

• an  increase  in  the  cathode  size, 

where  in  each  case  the  remaining  parameters  are  held  constant. 
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Table  3a 

Comparison  of  the  Total  Anodic  Current  Determined  by  Mathematical  and  Graphical 
Integration  for  a Fixed  Ratio  of  Anode  Radius  to  Cathode  Radius  a/c 


|q3  * to  ted , anodic 

4iroc(Ea  - Ec) 

Mathematical  Integration  Graphical  Integration 


0.0005 

0.252 

0.250 

0.318 

0.317 

0.375 

0.379 

0.01 

0.445 

0.445 

bulk 

0.467 

0.468 

0.001 

1.23 

1.17 

3.92 

| 

3.71 

bulk 

0.467 

0.468 

4.59 

4.56 

38.8 

38.1 

Table  3b 

Comparison  of  the  Total  Anodic  Current  Determined  by  Mathematical  and  Graphical 
Integration  for  Variable  Cathode  Size  with  a = 1,  b - 0.01,  £ = 10 
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APPENDIX  A 


MATHEMATICAL  SOLUTION  FOR  THE  ELECTRODE  POTENTIAL  OF 
CIRCULAR  CELLS  UNDER  BULK  ELECTROLYTE 


In  cylindrical  coordinates,  Laplace’s  equation 

3 2P(x,  y,  z)  + d2P(x,  y,  z)  + 32P(x,  y,  z) 


becomes 


dx‘ 


3 2P(r, 
3 r2 


3y" 


3 ^ 


(Al) 


1 3 PKr,z)  _ 3 2P(r,z)  n 

+ ^ — + — ns — _ u > 

r 3r  3 z 


(A2) 


where  the  9 term  is  deleted  because  P is  not  6 dependent.  The  general  solution  to  Eq. 
(A2)  is  well  known  but  will  be  outlined  here  for  the  convenience  of  the  reader.  One 
method  of  solution  is  separation  of  variables,  in  which  solutions  of  the  form 


P(r,z)  = R(r)Z(z) 


are  sought.  Thus 


3 P_ 
3 r 

3 2P 
3r2 


= R’(r)Z(z), 
= R"(r)Z(z) , 


(A3) 


( A4a) 


(A4b) 


and 


3 2P 
3Z2 


75  - R's)Z(z) , 


( A4c) 


where  the  prime  and  dot  indicate  differentiation  with  respect  to  r and  z respectively. 
Using  Eqs.  (A4)  in  Eq.  (A2)  and  making  appropriate  rearranrjaments  gives 


R"(r)  1 R'(r) 


R(r) 


+ — 


r R(r) 


Z(z) 

Z{z) 


( A5) 
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The  left  side  is  now  a function  of  r only  and  the  right  side  of  2 only,  so  that  each  side 
of  Eq.  (A5)  can  be  set  equal  to  an  arbitrary  constant,  -X2.  That  is, 

. -X2,  (A 

R(r ) r R(r) 


= -X2. 


The  first  separated  equation  is  thus 


R"(r)  + — R'(r)  + X2fl(r)  = 0.  (A8) 

After  making  the  substitution  s = Xr,  which  defines  s,  and  performing  some  algebraic 
manipulations,  Eq.  (A8)  takes  the  form 

d2R  1 dR  ( P2\ 

1?  + 717  + (^)*  = 0 fr"0*’  (A9) 

which  is  Bessel’s  equation  of  order  p and  has  the  independent  solutions  [Al,  A2]  R = J„(s) 
and  R - Yp(s),  where  Jp(s)  is  the  Bessel  function*  of  the  first  kind  of  order  p and  Yp(s) 
is  Neumann’s  Bessel  function^  of  the  second  kind  of  order  p.  The  general  solution  to  Eq. 
(A9)  is  then 


/?(r)  = AJ0(Xr)  + BY0(Xr) , (A10) 

where  A and  B are  constants.  The  solution  to  the  second  separated  equation 

Z(z)  - X2Z(2)  = 0 (All) 


Z(z)  = C\e^  + C'2e' 


(A12) 


2 (-lf 


(.)  - - T 

p \ 2/  L-‘_  n!(p  + n 


t s2  s*  l\  s6  / 1 l\ 

Vn(s)  = Jn( s)  log  s + -5  - -s — 9(1  + — + -r — 5 — , 1 . 

0 0 22  22*42  \ 2/  22-42-62\  2 3/ 
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where  C\  and  C2  are  constants.  Use  of  Eqs.  (A10)  and  (A12)  in  Eq.  (A3)  gives  the 
general  solution  of  Laplace’s  equation: 

P(r,z)  = [AJ0(Ar)  + BY0(Ar)]  [C\e^  + C'2e'X*]  . (A13) 

The  constants  A,  B,  C\ , and  C'2  are  to  be  evaluated  from  the  specific  boundary 
conditions,  which  have  been  listed  previously: 


(A14) 

(A15) 

(A16) 


r=c 

lim  P(r,z)  < JW 


P(r,  z)  - £ 

Differentiation  of  Eq.  (A13)  gives 
3 P(r,  z 


3 P(r,z) 

„ 3 z 

ii 

o 

II 

E , 0 < r < a ; 

a’  ’ 


E , a < r < c . 

C * 


(A17) 


3r 


= [-A  A«/j  (Ar)  - BAY’j(Ar)]  [Cje^  + C'2e'Xz] 


(A18) 


According  to  the  boundary  condition  of  Eq.  (A14),  Pr(r,  z)  = 0 when  r = 0,  so  that 

-AXJ^O)  - BAT1(0)  = 0,  (A19) 

which  requires  that  B = 0 because  1^(0)  is  infinite  and  ■/1(0)  is  already  zero  [A2].  If 
we  write  AC\  = Cj,  and  ,4C'2  = C2,  then 

P(r,z)  = JQ(\r ) [C^  +C2e-Xi]  . (A20) 

According  to  the  boundary  condition  of  Eq.  (A15),  Pr(r,  z)  = 0 at  r = c,  so  that 


-AJjtAcHC^*  + C2e-Xz]  = 0, 


(A21) 


or 


J,(Ac)  = 0.  ( A22) 

Let  these  roots  of  Jx  which  give  JjDO  = 0 be  called  as  shown  in  Fig.  Al:  that  is, 

= Arc  , n = 1,  2,  3,  ...  , ( A23) 
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or 

\ = — ■ (A24 

c 

Thus 


P(r,z)  = JQ  (xn  [CleXnZlc  + C2eXnZ/c]  . (A25 

The  boundary  condition  of  Eq.  (A16),  which  requires  P(r,  z)  to  be  bounded,  is  satisfied 
if  C j = 0.  Thus 


P(r,  z)  = C2</. 


0 («.  f) ' 


is  a particular  solution  to  the  problem.  The  general  solution  is  then  a linear  combination 
of  all  the  n solutions: 


oo 

P(r,  2)  = CQ  + T CnJ0(\r)e-*nZ, 


( A27 
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with  X„  = xnjc.  The  constants  Cn  are  evaluated  from  the  remaining  boundary  condition: 


P(r,  0)  - £ 


'dP(r,  z) 


dz 


= S(r)  = 


2=0 


£,  0 < r < a ; 


\Ec,  a < r < c . 


From  Eq.  (A27), 


P(r,  0)  = CQ  + £ CnJ0(\nr)  , 


n=l 


and 


- £ 


dP(r , *) 


3z 


^Ecnv0(V)- 

2=0  n=l 


Using  Eqs.  (A28)  and  (A29)  in  Eq.  (A17)  gives 

oo  OO 

C0  + £ C„W>  + £ *W0<V>  = s<r>’ 


n=l 


n=l 


or 


co  + £ c«(1  + f Xn>W>  = ^)- 


n = l 


One  next  employs  the  orthogonal  property  [ A3  ] : 

0,  n =£  m ; 


f r‘/0(Xnr)'/0<Xmr>dr  = | c2 

\ T [</°( 


(A17) 


(A28) 


(A29) 


( A30) 


( A31) 


( A32) 


Multiplying  both  sides  of  Eq.  (A31)  by  r</0(Xmr)  dr  and  integrating  from  r = 0 to  r = c 
gives  the  left-hand  side  (L.H.S.)  of  the  equation  to  be 

L.H.S.  = f C0rJo{\mr)dr  * f £ C„(l  + f dr  (A33) 

>=0  ->=0/1=1 

and  the  nght-hand  side  (R.H.S.)  to  be 


30 


NRL  REPORT  7835 


R.H.S.  = f S(r)rJ0(\mr)dr . 


The  first  integral  in  Eq.  ( A33),  denoted  as  /j : 


ll  = f Wa(kmr)dr, 


(A34) 


(A35) 


is  evaluated  first.  It  can  be  shown  [A2]  that 


JrJ0(Xmr)dr  = [rJ^r)]  , 


so  that 


but  J1(Xmc)  = 0 according  to  Eq.  (A22).  Thus  lx  = 0,  so  that 

L-H.S.  = f;  k(l  + £\)  f rJ0(\r)J^\mr)  dr  . 
n»l  L 0 J 

All  terms  are  zero  for  n ^ m,  so  the  only  term  which  survives  the  integration  is 

L.H.S.  - C„(l  ♦*>„)£  IW>12 
Operating  on  the  right-hand  side, 


(A36) 


(A37) 


( A38) 


( A39) 


R.H.S.  = J"  SirW^^r)  dr  + J S(r)rJQ(\m  r)  dr  , (A40) 


R.H.S.  = E_ 


r* 


rJ0^mr)dr  + Ec  f rJ0  <Xmr>dr- 
0 *r=a 


Use  of  Eq.  (A36)  gives 


R.H.S.  -is  [ftf,(Xma)l  + it  [cJ1(Xmc)-o^1(Xma)]  , 


(A41) 


(A42) 
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1 

but  again  J1(Xmc)  = 0;  thus  j 

’i 

R-H.S.  = a)  (£-£).  ( A43) 

m 

Equating  Eqs.  (A39)  and  (A43)  gives 

+ *\n>T  [■'ofVH2  = rJl  (K*HEa-Ec)-  <A44> 

m 

Solving  for  Cm  and  replacing  the  dummy  variable  m by  n gives  the  coefficients  Cn : 

J-J'{XnaUE°-EC) 

Cn  = — . ( A45)  ! 

One  next  solves  for  C0.  This  time  each  side  of  Eq.  (A31)  is  multiplied  by  r dr  and 
integrated  from  r = 0 to  r = c: 

c c °°  c • 

C0  f rdr  + f T C„(l  + nn)rJQ(\nr)dr  = [ S(r)  dr  . (A46) 

*V=  0 Jr=0  0 

Here  the  second  term  vanishes,  as  may  be  seen  by  evaluating  the  integral  according  to 
Eq.  (A36).  Again  the  right-hand  side  of  Eq.  (A46)  must  be  split  into  two  integrals.  The 
result  of  the  simple  integrations  in  Eq.  (A46)  is 


( A47) 


Use  of  the  expressions  for  Cn  and  C0,  as  given  by  Eqs.  (A45)  and  (A47),  in  Eq.  (A27) 
with  A„  = xn/c  gives  the  desired  expression  for  the  potential  distribution  in  the  system: 


P(r.z) 


E + (E  - E ) Y 

c a c ' L i 
n=  1 


( A48) 
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APPENDIX  B 

COMPUTER  PROGRAM  FOR  EVALUATION  OF  LOCAL  CURRENT  DENSITY  FOR 
COPLANAR  CONCENTRIC  CIRCULAR  ELECTRODES  UNDER  BULK  ELECTROLYTE 


PROGRAM  BJLK1 

c this  program  computes  the  local  current  density  along  the  radius 

c FOR  CO-PLANAR  CIRCULAR  ELECTRODES  COVERED  0Y  BULK  ELECTROLYTf 

c FOR  The  CASE  WHERE  ANODIC  ANO  CaThODIC  WAGNER  POLARIZATION 
C PARAMETERS  ARE  EQUAL. 

C 

c A»RADIUS  OF  ANODE 

C C-RADIUS  OF  CATHODE 

c R«DI STANCE  ALONG  RAOIJS 

C L«WAGNER  OOLARIZATION  PARAMETER 

C A/C«AR 

C L/C«LR 

C R/C*RR 

c x(ni«nth  zero  of  bessel  function  of  order  jord 

c 

REAL  LR.NJM2 

DIMENSION  X(5000) «F AC  TOR (5000) .TERM  (5000) 

•<•500 

c k-number  of  terms  in  sum 

READ  10.AR.LR 

10  FORMAT  (2X10*0) 

PRINT  11.AR.LR  „ ,,,, 

11  FORMAT  ( 1 Ml .5X« 7HAR«A/C*  Fl0.5.5*.7HLR»LZC*  F10.S.///I 

101  FORMAT^?25X.*U  ( x (N) 1 •»5X»8hX (N) *A/C. 7X.» J ( P) *. *X. 1 1 H ( A/C ) • J <P>> 
PRINT  102 

102  FORMAT  (26X,«0««27X»«1»»UX.«1*»>  ) 

PRINT  i o 3 

103  FORMAT  (6x.*N*»6X*#X(N)*i 9X . *8£SSEL2* . 9X  » »P»  » 9X » #0ESSEL 1 • » BX  . 
•*NUM2*,8X.»0EN0M1».8X» •0En0m2*. 7X . •FACTOR(N)*./ //) 

JORD* l 
NO*K 

CALL  8ESZER0( JORO.NO.X) 

DO  20  N«1.K 

P»X  < H ) #AH 

T«X (N)#LR 

8ESSEL1*0ESJ(P»i1 

num2»aR»0ESSEL1 

D£NOM1»1.0  ♦ T 

0ESSEL2»BESJ(X(N) .0) 

0£NOM2«8ESSEL2»»2 

FACTOR(N) «NUM2/ (D£NOM1«DENOm2) 

PRINT  iq.s.X(N) ,RESSEl2«P. BE SSELl.NUM2,DENOMl,OENOM2,F ACTOR (N) 

19  FORMAT  ( -X.I*.A(3X«F10.i»)  ,*(1X,E13.5)  ) 

20  CONTINUE 

READ  201 .RR.OELTA.RCUT 
201  FORMAT  ( F10.0) 

21  PRINT  22. RR 

22  FORMAT  <///3X«»RR»  »F5.3///) 

PRINT  221 
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221  FORMAT  (25X»8HX(N)«R/C*7X,«j  (*)•> 

PRINT  2->2 

222  FORMAT  (MX. •<>••/) 

PRINT  2?3 

223  FORMAT  (6X,»Na.7X,*X(N) *, i1X,»W*»9X,*BESSEL3#«6X.»FACT0R<N) *,BX, 
••TeHM(N)*.11X.*SUM«.///) 

SUM-0. 0 
00  30  N- 1 » K 
M»X(N1*RR 
BESSEL3-BESJ(N*0) 

TERM(N)  aF  ACTOR  (N)  ‘BESSEL3 
SUM-SUM. TERM(N) 

PRINT  29,N,X(N) ,w,8ESSEL3,FaCT0R(N) .TeRM(N) .Sum 

29  FORMAT  <3X,I*,1X,E12.5*2<3X,F10.5> . 3 ( 3X  ,E  1 . 5) ) 

30  CONTINUE 
RRaHH.OELTA 

IF  (RR.LE.RCUT)  GO  TO  21 
99  END 

SUBROUTINE  BES2£R0( JORO.NOt ZERO) 

IOENT  NUMBER  - C3007R08 

title  • zeros  of  the  bessel  function  of  the  first  kind 

IDENT  name  - C3-NRL 'BESZERO 
LANGUAGE  - 36002380  F8RTRAN 
COMPUTER  - CDC-3800 

CONTRIBUTOR  - JANET  P.  MASON.  CODE  7813 

RESEARCH  COMPUTATION  CENTER.  MIS  DIVISION 
ORGANIZATION  - NRL  - naval  RESEARCH  LABORATORY, 

KASHINGTOS,  n.C.  20390 
DATE  - 1 JULY  1971 

PURPOSE  - TO  FINO  Th£  FIRST  M ZEROS  OF  JSUBN(X)  FOR  o<N<5.  RhERE 
m IS  SUPPLIED*  pV  The  USER*  IN  The  SUBROUTINE  CALL 
DIMENSION  XJAOU) ,XJA1 (3) ,XJA2 (2) *XJA3(6) .XJA* (6) ,XJA5(9) .ZERO!  1) 
DA TA(XJA0a2.A0A8255577 ,5.520078 11 03* 8.6537279 129, 11.791 53**39) , 

1 (XJA1-3. 8317059702, 7. 0155866698, lo.l73*6Mi35)  , 

2 (XJA2.5. 1356223. 8. *172*41) , 

3 ( XJA3-6. 380 16 19 ,9, 7610231. 13.0152007.16.223*6*0, 19. *09*1*8. 

* 22.5827295), 

5 ( X JA*. 7. 5883*27, 1 1. 06*7o95,l*. 37  25367, 17. 6 159b60, 20. 8269330. 

6 2*. 0190195). 

7 ( XJA5. 8,77 1*838, 12.33860*2,15.70017*1 , 18.9801339,22.2177999, 

8 25,*3 03*1 1*28.6266183,31 .81 17167, 3*. 988781 3) 
pIa3. 1*15926536 

HOLOa*.o*JORO»JORO 

GO  TO  ( 1 ,2, 3, *,5,67 J0R0»1 

1 DO  11  1-1 , NO 

IF (I.GT.*)GO  TO  20 
ZEROlI)aXJAO(I) 

11  CONTINUE 
return 

2 DO  12  1-1, NO 

IF ( I , GT , 3 ) GO  TO  20 
ZERO ( I ) aX JA1 (I) 

12  CONTINUE 
RETURN 

3 00  13  1-1, NO 

IF ( I • GT , 2 ) GO  TO  20 
ZERO  < I ) aX  JA2 ( I ) 

13  CONTINUE 
return 

* DO  1*  1-1, NO 

IF ( I ,GT»6) GO  TO  20 
ZERO ( I ) aX  JA3 ( I ) 

1*  CONTINUE 
RETURN 


00000100 

0000^101 

00000102 

00000103 

0000010* 

00000105 

00000106 

00000107 

00000108 

00000109 

00000110 

00000111 

0000O112 

00000200 

00000300 

00000*00 

00000500 

00000600 

00000700 

00000800 

00000900 

00001000 

00001100 

00001200 

0000)300 

00001*00 

00001500 

00001600 

00001700 

00001800 

00001900 

00002000 

00002100 

00002200 

00002300 

00002*00 

00002500 

00002600 

00002700 

00002800 

00002900 

00003000 

00003100 

00003200 

00003300 

00003*00 
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5 00  15  I a 1 , NO 

IF ( 1 ,GT, 6) GO  TO  20 
ZERO t 1 ) aX  JA4 ( 1 ) 

15  CONTINUE 
RETURN 

6 DO  16  I ■ 1 • NO 
IF(I.GT.9)G0  TO  20 
ZERO ( I ) aX JA5 ( I ) 

16  CONTINUE 
RETURN 

20  BET*» (P 1/6.0) •<2.0#JO«O»4.0»I-l.Q> 

■ 1 »BET4»8 . 0 

w2bw1*m1 

ZERO ( I )■ BET*- ( HOLD* 1. 0> /» 1« ( 1.0*1. 0/*2»<4. 0*17. 0*HOLD-31.0) /3.0 

1 ♦1.0/«2#<32.0#  <83»0#H°LO*MOLO-982.0*ROLD.3779.0) /15.0 

2 *1 .0/*2» (66.0* ( 69*9. o*hOLO»HOLD«hOLD»1 53855, 0«HOLD*HOLD 

3 ♦ 15857*3. flaHOLO. 6277237.0) / 105.0) ) ) ) 

GO  TO  (11, 12. 13. 14,15,16) J0RD»1 

10  END 

FUNCTION  3ESJ ( X ,N) 

OATA(Ro«.28278**9»7E8) , (Hi»-.685265989lE7) , (R2-.38831312263E6) . 

1 (R3»-.9o57867*277E4) . (R4», 1Q8306963E3) , (R5a-.73  8.  335935) . 
2(R6«.29212672*87E-2) , < r.7a_.  65050  1 7fl57lE-5) , (R8a,6453801B051E-B) , 
3(So».28279449*7E8) , ( S 1 • . 2 1 695247743E6 ) , (S2a.7oO*6825l*7E3> , 

4 (Ao«2.5323*20902E2l . ( *1 ■* . 22 1 7704 1 l8El ) , ( A?a5 .24433 14672E-1 > , 

5 (B0-.44884594896E3I « ( 8 1 ■ . 75 322048S79E2 ) , 

6 (CO— 1. 233944555  IE  1) . <Cl»-2. 7 78892 1059) , (C2  — 4.9517399126E-2) , 
7(01».4100554523E2)  «(t,«64.)  . (G«4 .7  22366-.BE21 ) t 
8I00-.17496878239E3)  , 

A (RRO-.980  87274959E7) . ( HR  1 a-. 1 1 42532572 1 E 7 ) . ( HR2a . 4094621 3625E5) . 
9<RR3— .66660119856E3) , ( KH4* . 575754 1 4035E 1 ) , ( RR5»- . 279044755 1 9E- 1 
C(RR6». 734931321 11E*4)  , ( RB7-. 8430682 1 64 IE- 7)  , 

DISS0-.19617454991E8) , ( 851 « . 1 67 \ 1 673 1 84E6 ) , ( SS2».60777258247E3 > . 
E(0BO-.62836856631E3) . ( BB 1 « . 97 300 09462BE2 ) . 

F ( 000* . 2 1 1 9547833 1 E3 ) , (0U1..46917127629E2) . 

GIAA0-3.5451899975E2) , ( AA 1 .5 . 554434302  IE  1 ) » ( AA2»6.5223084285E-1 > , 
m(CC0«4.4822348228E1) t ( CC1 .9. 7 348068764 ) , ( CC2» 1 . 77255791 45E-1  > 
D»X*X 

IF(N.EQ.O)  GO  TO  6 * IF(N.EQ.l)  Go  TO  7 S GO  TO  B 

6 IF (D-F)  1,1,2 

1 Pa ( ( ( (RR*0»R7) «0*R6) *D*R5> *0»R4) #D  » P» ( ( (P»M3) »0»R2) »D*R1 ) »O»R0 
BESJ  »P/( ( (0»S2)«0*S1)«D*S0)  * RETURN 

2 IF(U.GT.G)  GO  TO  9 
A«ABS(X)  t 0“F/0 

Pa< (A2»0*Al)»0*A0)/( < 0*  B 1 ) *0  * B 0 ) 

Qa  ( (C2»3*C1)*D*C0)/(A*( (D.Ol ) •O.DO) ) 

BESJ  ■(COS(A)*(P*0>*SIN(A)«(P-Q) ) /SORT (A) 

RETURN 

7 IF (D«F >11,11,21 

11  Pa  ( < < < < (R97»0*RR6) *0*R85) *0.RR4) »0*RR3) •D»RRZ) #0*RR1 ) *D*RR0 
BESJ«X*P/( ( (0*SS2)*0*SS1)*0.5SO)  * RETURN 

21  IF(U.GT.G)  GO  TO  9 
A«A8S(X)  t 0-F/0 

Pa< (AA2»0»AA1) *0*AA0) /( (B.BBl) *0»BB0) 

Qa ( (CC2«0.CC1)»0«CC0) /(A»( (0*001 ) *D  *0D0 ) ) 

Aa ( COS(A) a(Q-P) »SIN(A) a(Q»P) ) /SQRT(A)  J I F ( X . L T . 0 > A»- A 
BESJaA 

return 

8 PRINT  81, N 

81  F0RMAT(//15X«E7R0R  IN  BESj,  N a#l5) 

GO  TO  100 

9 PR  I NT9 1 , X 

91  F0RHAT(//15X4ERR0R  IN  BE-Jt  ARGUMENT  X TOO  LARGE.  X a #E17.1o> 
100  BESJ«1 .E300 
EnO 
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APPENDIX  C 

COMPUTER  PROGRAM  FOR  EVALUATION  OF  LOCAL  CURRENT  DENSITY 
FOR  COPLANAR  CONCENTRIC  CIRCULAR  ELECTRODES 
UNDER  THIN  LAYERS  OF  ELECTROLYTE 


PROGRAM  THIN1 

THIS  program  computes  the  local  current  density  along  the  radius  I 

FOR  CO-PLaNaR.  concentric.  CIRCULAR  electrodes  COVERED  by  j 

A THIN  LAYER  OF  ELECTRBLYTE.  4 

i 

A»RADIUS  OF  ANODE 
C«RAOlUS  OF  CATHODE 
B«THICKNESS  OF  ELECTROLYTE  layer 
L-MAGNER  POLARIZATION  PARAMETER 
R-DISTANCE  ALONG  radius 
A/C*AR 
B/C*8R 
L/C*LR 
p/rspp 

X(N)«NTh  zero  OF  BESSEL  FUNCTION  OF  ORDER  J0«D 
anodic  and  cathodic  waoner  polarization  parameters  are  equal. 

real  LR.NJM1.NUM2 

DIMENSION  X (5000) • FACTOR (5000) * TERM(5000) 

K*500 

[ C K-Ni'MflER  OF  TERMS  IN  SUM 

I READ  I0.AR.BR.LR 

10  FORMAT  (3r10.0) 

PRINT  11,  AR,  BR,  LR 

11  FORMAT  U-U.5X.7HAR-A/C-  Fl0.5*5X»7HBR«B/C»  F10.5,5X,7HLR»L/C* 

•F10. 5.///) 

PRINT  lol 

101  FORMAT  (26X.PJ  (XN)**6X»6hXn*A/C*8X»*J  ( P) •» 8X, *TANH ( P) •,*X , 

••(A/C1J  (»)•) 

PRINT  102 

102  format  (2bx,»o*.25x»*i*»28x,*i»./> 

’ PRINT  103 

103  FORMAT  <6X.*N*,6X.*X(N)*i9X,*BESSEL2*.9X»»P*»9X,*BESSEL1**8X, 
••NUMl*.eX**NUM2*»8X,#DENDMl**8X»*DEN0M2**7X»*FA-T0R(N) •»///) 

J0RD«1 
NO-K 

CALL  BESZERO ( JORO. NO*  X) 

DO  20  Nal.K 
P«X(N)*AR 
C =X(N)*bR 
T«X(N)*i.R 

SINHQ«(eXP(Q)-EXP(»Q) 5/2.0 
COSHQ»(EXP(Q) *EXP(»Q) 5/2.0 


i 


'i 
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T~R~TTP^rrriwi|jpiinnT^iirpyiipiiiiiiwi^wpi  i n 


TaNHQ»SINh(«/COShQ 

NUMl«TANMa 

DENOH1-1.0  ♦ <T*7ANHC1) 

BESSEL  1 aBES J ( P» 1 ) 

BESSEL2aBESJ(X(N) ,0) 

NUM2»AR*BESSEL1 
0EN0M2* BESS EL 2**2 

FACTOR (N) ■ <NUM1*NUM2) / < DEnOmI *DEN0M2 > 

PRINT  1 9 , N , X ( N ) t BESSEL2tSf BESSEL l»NUMi ,NUM2, DENOMI , DENOM2.F ACTOR (N 

• ) 

19  FORMAT  (3X.I4.1X,E12.5.4(3X,F10.5) ,A(1X.E13.5, ) 

20  CONTINUE 

READ  201*RR. DELTA, RCUT 
201  FORMAT  (3F10.0) 

21  PRINT  22, RR 

22  FORMAT  (///3X.*RRa  *F5,3///) 

PRINT  221 

221  FORMAT  (25X.8HX (N) *R/C,7X,*J  <*)•> 

PRINT  222 

222  FORMAT  <*1X,*0*,/) 

PRINT  223 

223  FORMAT  (6X.*N*.7X.*X  <N)*» 1 1X,*W**9X.*BESSEL3*,6X.»FACT0R (N) *,8X. 
••TERM(N) «, 11X,*SUM*»///J 

SUM»0 • 0 
DO  30  Na 1 , K 
W»X  t N ) *RR 
BESSEL3aBESJ(W.O) 

TERM(N) aF ACTOR (N) *BESSEL3 
SUm»SUM*TeRM(N) 

print  29,n,x (n» ,w,bess6lj, Factor (n> , term < n> , sum 

29  FORMAT  < <X,I4,1X,E12.5,2(3X,F10.5) ,3<3X.E1 '.S> > 

30  CONTINUE 
RRaRR  ♦ DELTA 

IF  (RR.LE.RCUT)  GO  TO  21 
99  ENO 


SUBROUTINE  8ESZER0( JORO.NO, ZERO)  AND  FUNCTION  BESj(X,N)  aRE  GIVEN 
IN  APPENOIX  B. 
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